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Abstract 

We predict the branching ratios of Z — > e^/i , Z — ► e^r^ and Z — > /i^r 1 * 1 decays in the 
framework of the 2HDM with the inclusion of one and two extra dimensions, by considering 
that the new Higgs doublet is localized in the extra dimension with a Gaussian profile. 
We observe that their BRs are at the order of the mag nitude of 1(T 10 , 1(T 8 and 1(T 5 
with the inclusion of a single extra dimension, in the given range of the free parameters. 
These numerical values are slightly suppressed in the case that the localization points of 
new Higgs scalars are different than origin. 



*E-mail address: eiltan@heraklit.physics.metu.edu.tr 



1 Introduction 



The lepton flavor violating (LFV) interactions are sensitive tO the physics beyond the standard 
model (SM) and they are rich theoretically since they exist at least in the one loop level. The 
Z decays with different lepton flavor outputs, such as Z — > e/i, Z — > er and Z — > //r are among 
the candidates of LFV decays and they are clean in the sense that they are free from the long 
distance effects. In the literature, there is an extensive work on these decays PQ-JI|. The 
theoretical studies on such Z decays have been stimulated by the Giga-Z option of the Tesla 
project which aims to increase the production of Z bosons at resonance. 

In the framework of the SM the lepton flavor is conserved and its extension with massive 
neutrinos, so called i/SM model, permits the LFV interactions with the lepton mixing mecha- 
nism ^3]. However, in this model, the theoretical predictions of the branching ratios (BRs) of 
these LFV Z decays are extremely small when the internal light neutrinos are light PJ Ej 

BR(Z -> e ± p ± ) ~ BR(Z -> e ± r ± ) ~ 10" 54 , 

BR(Z -> jjV 1 ) < 4 x 1(T 60 . (1) 

These numbers are far from the experimental limits obtained at LEP 1 [3]: 
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and from the improved ones at Giga-Z [7j: 



BR(Z -> e^) < 2 x 1(T 9 , 
BR(Z -> e ± r ± ) < / x 6.5 x 1(T 8 , 

^^(Z -> // ± r ± ) < / x 2.2 x (3) 
with / = 0.2 — 1.0. Here the BRs are obtained for the decays Z — > Z1Z2 + hh, namely, 

BR(Z _ = r(Z ^^ 2 + [2/l) . (4) 

t z 

The extensions of z/SM with one heavy ordinary Dirac neutrino |2], and with two heavy 
right-handed singlet Majorana neutrinos [2] ensure to enhance the BRs of the corresponding 
LFV Z decays. The possible enhancements in their BRs have been obtained in other models; 
in the the Zee model jHj, the two Higgs doublet model (2HDM), without (with) the inclusion 



of the extra dimension [HK^U]), the supersymmetric models [IIJEI| an d the top-color assisted 
technicolor model [T3j . 

In this work, we study the LFV processes Z — > e ± /x ± , Z — > e ± r ± and Z — > fi r in the 
framework of the 2HDM with the inclusion of a single (two) extra dimension(s). Here the LFV 
interactions are induced by the internal new neutral Higgs bosons h° and A at least in the 
one loop level. The extension of the Higgs sector brings new contribution to the BRs of the 
considered decays. On the other hand, the inclusion the extra dimensions enhance the BRs 
since the particle spectrum extends after the compactification of the extra dimensions. 

The extra dimension idea was originated from the study of Kaluza-Klein [16J which was 
related to the unification of electromagnetism and the gravity and the motivation increased 
with the study on the string theory which was formulated in a space-time of more than four 
dimensions. Since the extra dimensions are hidden to the experiments at present, the most 
favorable description is the compactification these new dimensions to the surfaces with small 
radii. In the case of that the extra dimensions are at the order of submilimeter distance, for two 
extra dimensions, the hierarchy problem in the fundamental scales could be solved and the true 
scale of quantum gravity would be no more the Planck scale but in the order of electroweak 
(EW) scale [IZl HE! - 

The effects of extra dimensions on various physical processes have been studied in the 
literature extensively JTj-jlB]- ^ n ^ ne ex ^ ra dimension scenarios, the compactification procedure 
leads to appear new particles, namely Kaluza-Klein (KK) modes in the theory. If all the fields 
feel the extra dimensions, so called universal extra dimensions (UED), the extra dimensional 
momentum, therefore the KK number at each vertex, is conserved. If some fields feel the extra 
dimensions but not all in the theory, those extra dimensions are called non-universal extra 
dimensions and this is the case where the KK number at vertices is not conserved. The non- 
conservation of the KK number at the vertex results in the possibility of the existence of the 
tree level interaction of KK modes with the ordinary particles. In another scenario, some fields 
are considered to be localized in the extra dimension(s). In the split fermion scenario pl2~]-[39j. 
the fermions are assumed at different points in the extra dimension with Gaussian profiles and 
this ensures a possible solution to the hierarchy of fermion masses by considering the overlaps 
of fermion wave functions in the extra dimensions. The localization of the Higgs doublet in 
the extra dimension has been considered in [40J, by introducing an additional localizer field. 
In |4T] the branching ratios of the radiative LFV decays have been studied in the split fermion 
scenario, with the assumption that the new Higgs doublet is restricted to the 4D brane or to 
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a part of the bulk in one and two extra dimensions, in the framework of the 2HDM. 02] is 
devoted to analysis of the BRs of the radiative LFV decays in the case that the new Higgs 
scalars were localized in the extra dimension with the help of the localizer field and the SM 
Higgs was considered to have a constant profile. In the recent work [13], the radiative LFV 
decays were studied with the assumption that the new Higgs doublet was localized in the extra 
dimension with a Gaussian profile, by an unknown mechanism, however, the other particle zero 
modes have uniform profile in the extra dimension. 

The present work is devoted to the BRs of the LFV Z decays in the 2HDM, with the inclusion 
of one and two extra dimensions by considering that the new Higgs doublet is localized in the 
extra dimension with a Gaussian profile, by an unknown mechanism, however, the other particle 
zero modes have uniform profile in the extra dimension. First, we assume that the new Higgs 
doublet is localized around origin and, second, we take the localization point as different than 
the origin but near to that. We observe that the BRs of the LFV Z decays Z — * e ± /i ± , 
Z — ► e ± r ± and Z — ► /i ± r ± reach to the values at the order of the magnitude of 10~ 10 , 10~ 8 and 
10 -5 with the inclusion of a single extra dimension, in the given range of the free parameters. 
These numerical values are slightly suppressed in the case that the localization points of new 
Higgs scalars are different than origin. 

The paper is organized as follows: In Section 2, we present the effective vertex and the BRs 
of LFV Z decays in the 2HDM with the inclusion of extra dimensions. Section 3 is devoted 
to discussion and our conclusions. In appendix section, we give the explicit expressions of the 
factors appearing in the effective vertex. 

2 The effect of the localization of the new Higgs doublet 
on the lepton flavor violating Z —> Ifl^ decays in the 
framework of the two Higgs Doublet model. 

The LFV Z boson decays Z — ► exist at least in the one loop level and, therefore, the 
theoretical values of the BRs are extremely small in the SM. With the extension of the Higgs 
sector in which the flavor changing neutral current (FCNC) at tree level is permitted, there 
appear additional contributions to the BRs of the LFV processes. The multi Higgs doublet 
models are among the candidates for such models and, in the present work, we consider the 
2HDM with FCNC at tree level. The inclusion of the the spatial extra dimension further 
enhances the BRs, since the particle spectrum is extended after its compactification and the KK 
modes of the fields which are accessible to the extra dimension bring additional contributions. 
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Here, the idea is to consider that the new Higgs scalars are localized in the extra dimension, 
with Gaussian profiles, by an unknown mechanism, and, the other particles have constant zero 
mode profiles in the extra dimension. 

The Yukawa Lagrangian responsible for the LFV interactions in a single extra dimension 
reads, 



C Y = tfJiLfaEjR + h.c. , 



(5) 



where L and R denote chiral projections, L(R) = 1/2(1 ^75), 02 is the new scalar doublet and 
are the FV Yukawa couplings in five dimensions, where i,j are family indices of leptons, k 
and Ej are lepton doublets and singlets respectively. These fields are the functions of and 
y, where y is the coordinate represents the fifth dimension. Here we choose the Higgs doublets 

0x and 02 as 
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and their vacuum expectation values as 
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This choice makes it possible to collect the SM (new) particles in the first (second) doublet 
and Hi and H 2 becomes the mass eigenstates h° and A , respectively since no mixing occurs 
between two CP-even neutral bosons H° and h° at tree level. 

The five dimensional lepton doublets and singlets and the SM Higgs field are expanded into 
their KK modes with the compactification of the extra dimension on an orbifold S 1 /Z 2 with 
radius R and they read 



0i (x, y) = 
k(x,y) = 
Ei(x,y) = 



1 



V2^R 
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x) + V2Y / (p i r\x)cos(ny/R)\ , 

n=l J 

00 

if^x) + V2J2 [&(x) cos(ny/R) + sm(ny/R) 

n=l J 

00 

E^(x) + y/2 E [E&\x) cos(ny/R) + E$\x) sin(ny/R) 



(8) 



where (^^(x), lfi(x) and E^(x) are the four dimensional Higgs doublet, lepton doublets and 
lepton singlets respectively. On the other hand the new Higgs scalar profiles read, 



S(x,y) = Ae- py2 S(x), 



(9) 
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and the mechanism beyond the localization is unknown 1 . Here the normalization constant A is 

A= PA"* . (10) 

7rV4 JErfly/TfiiiR] 



The strength of the localization of the new Higgs doublet in the extra dimension is regulated 
by the parameter (3 = l/a 2 , where a, a = pR, is the Gaussian width of S(x,y) in the extra 
dimension. Here the function Erf[z] is the error function, which is defined as 

Erf[z] = e~ t2 dt . (11) 



"K JO 

The modified Yukawa interactions in four dimensions can be obtained by integrating the com- 
bination fi^(R){ x ,y) S(x,y) fj 7 R(L)(^,y), appearing in the Lagrangian (eq. J5J)), over the fifth 
dimension as 

I = *V f8!$frv) S(x,y) fj$$(x,y) , (12) 
where fj^nl are the KK basis (zero-mode) for lepton fields (eq. (JBJl), and finally, we get 

I = V n f$$(x)S(x)fji$(x), (13) 

with the factor V n 

V n = Ac n , (14) 
and the function A which is defined in eq. (JTUJ). The function c n in eq. (|T4*jl is obtained as: 



.,2 



c n = V t ^ • (15) 

Notice that the factor A is embedded into the definition of the the Yukawa couplings £^ in four 
dimensions as 

e5=A&v (i6) 

where are the Yukawa couplings in five dimensions (see eq. (J3j)). 

In the following, we consider that the new Higgs scalars are localized in the extra dimension 
at the point y H , y H = a a near to the origin, namely, 

S{x,y) = A H e-K*-*rf S(x) , (17) 



1 We consider the zero mode Higgs scalars and we do not take into account the possible KK modes of Higgs 
scalars since the mechanism for the localization is unknown and we expect that the those contributions are 
small due to their heavy masses. 
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with the normalization constant 



Ah = *M!1 (is) 

(2tt)V4 J Erf [y/2]3 (nR + y H )] + Erf[y/2]3 (nR-y H )]' 



The integration of the combination fjxPm(x,y) S(x,y) fj^ul(x,y) over extra dimension brings 
the factor V n appearing in eq. (fT3Jl as 

V n = A H c n , (19) 
where Ah is the normalization constant defined in eq. (jl8j) and function c n reads 



Similar to the previous case, we define the Yukawa couplings in four dimensions as 

^ = A H ^. (21) 

Now, we would like to make the same analysis in the case of two spatial extra dimensions. The 
six dimensional lepton doublets and singlets and the SM Higgs fields are expanded into their 
KK modes with the compactification of the extra dimension on an orbifold (S 1 x S 1 )/Z 2 , with 
radius R and they read 



l — Uf'°\x) + 2j2<P ( r\x) <xx(ny/R + 8z/R)} , 

k(x, y, z ) = -L U° L °\x) + 2 J2 \liF\z) cos(ny/R + sz/R) + l\ n R s \x) sm(ny/R + sz/R) 

Ei{x, V,z) = -^ &° R °\x) +2 £ [E&'\x) cos(ny/R + sz/R) +E$'\x) sm(ny/R + sz/R)] \ , 

(22) 

where ^'^(x), lfif\x) and Ef R \x) are the four dimensional Higgs doublet, lepton doublets 
and lepton singlets respectively. Here the summation is done over the indices n, s but both 
are not zero at the same time. Similar to a single extra dimension case, the new Higgs scalar 
profiles read, 

S(x, y, z) = A'e-^ y2+z2) S(x) , (23) 
and the mechanism beyond its localization is unknown. Here normalization constant A' is 

A' = {2 ^' 2 (24) 

irWErfiy/TPicR]' 1 ' 
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The modified Yukawa interactions in four dimensions can be obtained by integrating the com- 
bination fiL(R\ \x,y) S(x,y,z) fjR(L\ {x,y) over the fifth and sixth dimensions: 

I = dy r dz f$ffi\x, y, z) S(x, y, z) f$$f>\x, y, z) , (25) 

where 

I = Vn,s fiL^''^ (x) S(x) /j^'(|)' 0)) (x) , (26) 

with the factor V n>s 

Vn,s = A c n>s , (27) 
and the function A' which is defined in eq. (|24j) . The function c„ iS in eq. (}2Tj) is obtained as: 

= ~ mhr- -' (28) 

Here the Yukawa couplings £y in four dimensions read 

c/; (29) 

where £jf^- are the Yukawa couplings in six dimensions. 

The internal neutral Higgs particles h° and A play the main role in the existence of the 
Z — > l^l£ decay, theoretically In Fig. dthe necessary 1-loop diagrams, the self energy and 
vertex diagrams, are given. The inclusion of extra dimensions brings additional lepton KK 
mode contributions. The general effective vertex for the interaction of on-shell Z-boson with a 
fermionic current reads 

1^ = 7 M (/y - Ia 75) H Um + !e 75) cruu q u , (30) 

mw 

where q is the momentum transfer, q 2 = (p — p') 2 , fy (Ja) is vector (axial-vector) coupling, 
fhi (/b) magnetic (electric) transitions of unlike fermions. Here p (— p') is the four momentum 
vector of lepton (anti-lepton). The vector (axial- vector) fy (f a) couplings and the magnetic 
(electric) transitions fu (/e) including the contributions coming from a single extra dimension 
can be obtained as 

liV ) ' 
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/* = EOS^E/^), ( 31 ) 



i=l n=l 

where I^vame) are ^he couplings without lepton KK mode contributions and they can be 
calculated by taking n = in eq. On the other hand the couplings fl?v,A,M,E) are the 

ones due to the KK modes of the leptons (see eq. (f35J)). Here the summation over the index 
i represents the sum due to the internal lepton flavors, namely, e,/i, r. We present f^y a m e) 
in the appendix, by taking into account all the masses of internal leptons and external lepton 
(anti- lepton). If we consider two extra dimensions where all the particles are accessible, the 
couplings fi{yAME) appearing in eq. (J3T]) should be replaced by f^y\ m e) anc ^ they reac ^ 

3 oo 

fv = E(# 0) + 4E# 

i=l n,s 
3 oo 

f A = E(/S' 0) + 4E/S' 

i=l n,s 
3 oo 

Im = E(A 0) + 4E^ (n ' s) ; 

i=l n,s 
3 oo 

f E = E(/S 0) +4E/£' S) )> (32) 

i=l n,s 

where the summation would be done over n, s = 0, 1, 2... except n = s = (see appendix for 
their explicit forms). 

Finally, the BR for Z — > V[ l£ can be written in terms of the couplings fv, /a, /at and Je 

as 

BH(Z ir ij) = ^ + + ^(I/mI 2 + l/.l 2 )} (33) 

where is the total decay width of Z boson. In our numerical analysis we consider the BR 
due to the production of sum of charged states, namely 

BR{Z -> If If) = r ( Z ->fck + Mi) (34) 

3 Discussion 

The LFV Z decays Z — > Z]* 1 /^, Zi 7^ Z 2 are rare decays in the sense that they exist at least in 
the one loop level and they are rich theoretically since the physical parameters of these decays 
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contain number of free parameters of the model used. In the framework of the 2HDM the 
internal leptons and new scalar bosons drive the interaction and the corresponding physical 
quantities are sensitive to the Yukawa couplings 2 i, j = e,\i,T, which are among the free 
parameters of the model. These couplings should be restricted by using present and forthcoming 
experiments. Here, we assume that the couplings which contain r index are dominant respecting 
the Cheng-Sher scenario 0^ and, therefore, we consider only the internal r lepton in the loop 
diagrams. In addition to this, we take the Yukawa couplings as symmetric with respect 
to the indices % and j. As a result, among the Yukawa couplings we need the numerical values 
for the £jy >Te , Cntu anc ^ £n,tt- Furthermore, the new Higgs masses are also free parameters of 
the model and we take their numerical values as m^o = 100 GeV, ttiao = 200 GeV. 

In the present work, we study the LFV decays Z — > if if, l\ ^ h in the framework of 
the 2HDM with the addition extra dimensions. Our assumption is that the new Higgs scalars 
are localized in the extra dimension with Gaussian profiles by an unknown mechanism, how- 
ever, the other particles have uniform zero mode profiles in the extra dimension. Here we 
choose one (two) extra dimension(s) which are compactified on to orbifold S 1 /Z 2 ((S 1 x S 1 )/^) 
with the compactification scale 1/R, which is another free parameter. The direct limits from 
searching for KK gauge bosons imply 1/R > 800 GeV, the precision electro weak bounds 
on higher dimensional operators generated by KK exchange place a far more stringent limit 
1/R > 3.0 TeV [22] and, from B — > (j)Ks, the lower bounds for the scale 1/R have been ob- 
tained as 1/R> 1.0 TeV, from B — > ip K$ one got 1/R> 500 GeV, and from the upper limit 
of the BR, BR (B s -> < 2.6 x 10~ 6 , the estimated limit was 1/R> 800 GeV jHSJ- On 

the other hand, the localization of new Higgs doublet is regulated by the parameter a, which 
is the Gaussian width of new Higgs doublet in the extra dimension, and it is chosen so that 
it does not contradict with the experimental results. Here, we take the compactification scale 
1/R in the range 200 GeV < 1/R < 1000 GeV and choose the Gaussian width a = p R at most 
0.05 R. Notice that throughout our calculations we use the input values given in Table ((TJ). 



2 In the following we use the dimensionful coupling £jy i ■ in four dimensions, with the definition £ 




^~7n~ £n a wnere N denotes the word "neutral" . 



9 



Parameter 


Value 


rrifi 
m T 
m w 
m z 

m^o 

G F 

Tz 

sin 9 W 


0.106 (GeV) 
1.78 (GeV) 
80.26 (GeV) 
91.19 (GeV) 
inn ffipv^ 

1UU I vac V J 

200 (GeV) 
1.1663710- 5 (GeV- 2 ) 
2.490 (GeV) 
V0.2325 



Table 1: The values of the input parameters used in the numerical calculations. 

In our analysis, we first consider that the new Higgs doublet is localized around the origin in 
a single extra dimension. Furthermore, we choose the localization point is near to the origin, at 
the point yu = a a, and study its effect on the BRs. We continue to analyze the same physical 
quantity with the inclusion of two extra dimensions. 

Fig. El is devoted to the parameter p = a/R dependence of the BR (Z — > p ± e^) for 
£# iTe = 0.1 GeV, £% tTli = 10 GeV and l/R = 500 GeV. Here the solid (dashed, small dashed) 
line represents the BR without lepton KK modes in a single extra dimension (with lepton KK 
modes in a single extra dimension, with lepton KK modes in two extra dimensions). It is 
observed that BR is at the order of the magnitude of 10 -14 without lepton KK modes in a 
single extra dimension, for the the parameter p ~ 0.05. In this case the BR is sensitive to p. 
With the inclusion of lepton KK modes, the BR enhances to the values of the order of 10~ 10 
and this is almost four order enhancement in the BR. For two extra dimensions, the numerical 
value of the BR is slightly smaller compared to the single extra dimension case, since there is an 
additional suppression factor (see the exponential factor in eq. (J2BJ)) appears in the expressions. 

In Fig. El we present the p dependence of the BR [Z — > r ± e^) for £^ Te = 0.1 GeV, 
£jv rr = 100 GeV and l/R = 500 GeV. Here the solid (dashed, small dashed) line represents 
the BR without lepton KK modes in a single extra dimension (with lepton KK modes in a 
single extra dimension, with lepton KK modes in two extra dimensions). This figure shows 
that BR is at the order of the magnitude of 10~ 12 without lepton KK modes in a single extra 
dimension, for the the parameter p ~ 0.05 and its sensitivity to the parameter p is strong. The 
inclusion of lepton KK modes results in an considerable enhancement in the BR almost four 
order and the BR of the decay under consideration enhances to the values of the order of 10~ 8 . 
For two extra dimensions, the numerical value of the BR is almost the same as the single extra 
dimension case. 
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Fig. |U represents the p dependence of the BR (Z — > p~^) for £jv T „ = 10 GeV, £j} TT = 
100 GeV and 1/R — 500 GeV . Here the solid (dashed, small dashed) line represents the BR 
without lepton KK modes in a single extra dimension (with lepton KK modes in a single extra 
dimension, with lepton KK modes in two extra dimensions). We observe that BR is at the order 
of the magnitude of 10~ 8 (10~ 5 ) without (with) lepton KK modes in a single extra dimension, 
for the the parameter p ~ 0.05. Similar to the previous decays, the inclusion of lepton KK 
modes results in an considerable enhancement in the BR more than three orders. For two extra 
dimensions, the numerical value of the BR is almost the same as the single extra dimension 
case. 

At this stage, we study the dependence of the BR of the LFV Z decays to the Yukawa 
couplings, regulating the lepton-lepton-new Higgs interactions. Figs. EHEK represent the 
Yukawa coupling £jf Te - Cntb ~ £n,tt dependence of the BR [Z — > /i ± e ± )-BR (Z — > p ± e ± )- 
BR [Z -> r^fj^) for g> T|1 = 10GeV-£% tTT = 100 GeV = 10 GeV, p = 0.01 and 
1/R — 500 GeV. Here the solid (dashed, small dashed) line represents the BR without lepton 
KK modes in a single extra dimension (with lepton KK modes in a single extra dimension, 
with lepton KK modes in two extra dimensions) for all figures. The BR is strongly sensi- 
tive to to the Yukawa coupling Te - £,§ Te - TT an d m the interval 0.005 < £§ Te < 0.05- 
0.005 <lfj Te < 0.05-50 <lff TT < 100 it enhances almost two-two -one order of magnitude. 
These figures also show that the inclusion of lepton KK modes causes the BR to increase 
considerably. 

Finally, we study the effects of the position of the localization point of the new Higgs doublet 
on the BR of the considered decays. 

Fig. IHlrepresents the parameter a dependence of BR [Z — > p^ 1 e ± ) for the Yukawa couplings 
= io Gel/, ff, re = 0.1 GeV, p = 0.01 and 1/R = 500 GeV. Here the solid (dashed) line 
represents the BR without lepton KK modes in a single extra dimension (with lepton KK 
modes in a single extra dimension). Without lepton KK modes, the BR is not sensitive to the 
parameter a for the interval 0.1 < a < 1. The inclusion of lepton KK modes makes the BR 
sensitive to the parameter a and increasing values of this parameter results in to decrease the 
BR almost one order for the interval taken. 

Fig. IUHTU1 represents the parameter a dependence of BR [Z — > r ± e^ - BR [Z r ± p^ 1 ) 
for the Yukawa couplings £f TT = 100 GeV, £f jTe = 0.1GeV-^ Tfl = 10 GeV, p = 0.01 and 
1/R — 500 GeV. Here the solid (dashed) line represents the BR without lepton KK modes in 
a single extra dimension (with lepton KK modes in a single extra dimension) for both figures. 
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The BR is not sensitive to a for the interval 0.1 < a < 1 without lepton KK modes for both 
decays. The inclusion of lepton KK modes increases the sensitivity of the BR to the parameter 
a and increasing values of this parameter results in to decrease the BR almost one order for 
the considered interval for both decays. 

As a summary, the BR (Z — > /i ± e ± ) ( (Z — > re), (Z — > r ± /i ± )) enhances up to the 
values of the order of 10~ 10 (10 -8 , 10 -5 ) with the inclusion of lepton KK modes in a single 
extra dimension. For two extra dimensions, the numerical value of the BRs are slightly smaller 
compared to the single extra dimension case. On the other hand the inclusion of lepton KK 
modes makes the BRs sensitive to the parameter a and increasing values of this parameter re- 
sults in to decrease the BR almost one order for the considered interval of this parameter. With 
the forthcoming more accurate experimental measurements of the these decays, the valuable 
information can be obtained to detect the effects due to the extra dimensions and the possible 
localization of the Higgs doublet. 

A The explicit expressions appearing in the text 

Here we present the explicit expressions for fjy\ fj$, fj$ and f^' (see eq. (JHU): 
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o r tier 9 r tier 



1 r uer r ver 

/„ -^AQ feO -kfeQ AQ 

gm w 



J q dx jf dy j ^(1 - x - y) (cy r/ y + c A r/f ) (x m,- + y m,+ ) 



64 7T 2 COS ^jy 
777,^ (C A (X - y) I]' + Cy 77+ (x + y)) 



T ver 



^(1 - x - y) (cvvY + °Avt) {xm^ +ym l +) -m, (c A (x - y) ^ +c v nf (x + y))^J 



(1-x-y) ( Vi {X ^" + V ^ (— + + ^ f 1 

V » / I o \ r «er r tier / o \ 



r ver j ver / o \ r ver T ver 



gm w rl rl ~ x 



64 7r 2 cos 9 W Jo 



J Q dx dyll(l-x-y)(-{c v nf + c A nY){xm l --ym l +) S ) 



mi (c A (x - y) 77+ + cy 77^ (x + y)) 



-k/iO 



^(1 -x -y) (- (c v nf + c A nY) (xm,- -ym,+ )) + mi (c A (x - y) r/+ + cy ^ (x + y))^J 
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+ (1-x-y) 

where 



(m-x-m + V ) (— + — — VU 

V ii IX <* ' \ Tver Tver I 9 V Tver Tver I I f ' 



L%h° = m h° (1 — X ) + ( m ! n)2 _ m f- (1 — X )) ^ ) 

r sel f T sel f I \ 

L i,ao = L x ^{m h a -> m A o) , 

T sel f T sel f i \ 

l 2,ho = L iM m ii ^ m i+)' 

j self j self i \ 

U 2, A° — ^l^OK 771 ^ ~^ m l+) ' 

tffl = rnla (1 - x - y) + mf' 2 (x +y) - q 2 xy , 

L v h t r A0 = m 2 A0 x + mf )2 (1 - x - y) + (m 2 h0 -q 2 x)y, 

L\% = L v h l r (m h o -> m A o) , 

LXo\o = L^ AO {m h o -> m^o) , (36) 



and 



Vi C n {^ih^il2 ^hi ' 

^ = C n {£zfi £if 2 + %>Ui£hi } ' 

% = c n {£zfi £if 2 — Cfi^fi} • (37) 

The parameters cv and c A are = — | and cy = | — sin 2 6 W and the masses mf^ read 
= \Jm 2 + n 2 /R 2 , where i? is the compactification radius. In eq. ()37|) the flavor changing 
couplings represent the effective interaction between the internal lepton i, (i = e, /i, r) 
and outgoing (incoming) j = 1 (j = 2) one. The parameter c n is defined in eq. (j!5|) for the 
localization of the new Higgs doublet around the origin and in eq. (j2*Uj) for the localization 
of the new Higgs doublet around the point y# near to the origin. In the case of two extra 
dimensions c n is replaced by c njS (see eq. (J2EJ)) and the masses m\ are replaced by m\ n ' s \ 



mf - ^ = \Jm 2 + m 2 n + m 2 , with m n = n/R, m s = s/R. 

Finally, the couplings £^ may be complex in general and they can be parametrized as 

e£- = l^|e^ , (38) 

where i, lj denote the lepton flavors and % are CP violating parameters which are the possible 
sources of the lepton EDM. However, in the present work we take these couplings real. 
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(c) (d) 



Figure 1: One loop diagrams contribute to Z — > k + j~ decay due to the neutral Higgs bosons 
ho and A in the 2HDM. % represents the internal, j (k) outgoing (incoming) lepton, dashed 
lines the vector field Z, h and A fields. In 5 (6) dimensions there exist also the KK modes of 
lepton and Higgs fields. 
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p 



Figure 2: BR( Z — > e^) with respect to the parameter p for £^ Te = 0.1 GeV, ^ Tfl = 10 GeV 
and 1/R = 500 GeV. Here the solid (dashed, small dashed) line represents the BR without 
lepton KK modes in a single extra dimension (with lepton KK modes in a single extra dimension, 
with lepton KK modes in two extra dimensions). 
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Figure 3: BR( Z -> r ± e ± ) with respect to the parameter p for £^ Te = 0.1 GeV, £^ rT = 100 (7eV 
and 1/R = 500 GeV. Here the solid (dashed, small dashed) line represents the BR without 
lepton KK modes in a single extra dimension (with lepton KK modes in a single extra dimension, 
with lepton KK modes in two extra dimensions). 
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Figure 4: BR( Z — >■ t± /i ± ) with respect to the parameter p for ^ Tfl = 10 GeV, £$ jTT = 100 GeV 
and 1/R = 500 GeK. Here the solid (dashed, small dashed) line represents the BR without 
lepton KK modes in a single extra dimension (with lepton KK modes in a single extra dimension, 
with lepton KK modes in two extra dimensions). 
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Figure 5: Z -> /x ± e ± with respect to f$ jTe for ^ T/i = 10 GW, p = 0.01 and = 500 GeV. 
Here the solid (dashed, small dashed) line represents the BR without lepton KK modes in a 
single extra dimension (with lepton KK modes in a single extra dimension, with lepton KK 
modes in two extra dimensions). 
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Figure 6: Z -> r ± e ± with respect to £f jTe for ^ TT = 100 GeV, p = 0.01 and 1/i? = 500 GeV. 
Here the solid (dashed, small dashed) line represents the BR without lepton KK modes in a 
single extra dimension (with lepton KK modes in a single extra dimension, with lepton KK 
modes in two extra dimensions). 
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Figure 7: Z -> r ± /i ± with respect to f$ |TT for = 10 Gel/, p = 0.01 and 1/R= 500 Gel/. 
Here the solid (dashed, small dashed) line represents the BR without lepton KK modes in a 
single extra dimension (with lepton KK modes in a single extra dimension, with lepton KK 
modes in two extra dimensions). 




0.001 



1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



a 



Figure 8: Z — > /x ± e ± with respect to the parameter a for £§ Tfl = 10 Gel/, £§ Te = O.lGeV, 
p = 0.01 and 1/R = 500 GeV. Here the solid (dashed) line represents the BR without lepton 
KK modes in a single extra dimension (with lepton KK modes in a single extra dimension). 
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Figure 9: Z — > e ± with respect to the parameter a for £§ TT = 100 GeV, £§ Te = 0.1 GeV, 
p = 0.01 and 1/R = 500 GeV. Here the solid (dashed) line represents the BR without lepton 
KK modes in a single extra dimension (with lepton KK modes in a single extra dimension). 
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Figure 10: Z — > fi^ 1 with respect to the parameter a for Cn,tt = 100 GeK, £§ T „ = 10 GeV, 
p = 0.01 and 1/R = 500 GeV. Here the solid (dashed) line represents the BR without lepton 
KK modes in a single extra dimension (with lepton KK modes in a single extra dimension). 
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